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The density operator in the Luttinger model consists of two components one of which describes 
long-wave fluctuations and the other, charge density wave (CDW) component, describes short- 
range electron correlations. It is commonly believed, that the conductance is determined by the 
smooth component, and the CDW component is addressed to only when an impurity is present. 
We investigate the contribution of the CDW component to the dynamic density response of a one- 
dimensional system free from any impurities. The conventional form of the CDW density operator is 
shown not to conserve the number of particles in the system. The corrected CDW density operator 
devoid of this shortcoming is proposed. The contribution of the CDW component to dissipative 
\Q • conductance is found to dominate over that of the long-wave component in the case when the Id 

conductor is locally disturbed by a low-frequency potential. 
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Electron transport in quantum wires is of significant interest because of the fundamental role that electron-electron 
interaction plays there. It is known that a|H.-arbitrarily small interaction destroys the quasi-particle excitations in a 
one-dimensional (Id) normal Fermi~systeml!m The excitations of an interacting Id system are of bosonic nature. They 
represent the waves of electron density, transferring charge and producing current. These waves have two components: 
q ■ one is smooth on the scale of the Fermi wavelength, and the other rapidly oscillates in space. 

The contribution of the smooth density component to transport is now well investigated.LlTj However, the theoretical 
predictions pertaining to this part of the density are not unambiguously supported by experiment to date. Moreover, 
there is a principal question of what effects, connected inherently with electron-electron interaction, are observable 
in transport. Turning to this question, we draw attention to the fact that the Luttinger liquid is a strongly corre- 
lated electron state. Therefore, exploring interaction effects inevitably must involve the rapidly oscillating density 
, component, which describes short-range electron correlations. 
OS ■ This component of electron density comes from the interference of electrons, moving in opposite directions. It 
resembles the charge density waves (CDW) in Peierls systems, but, in contrast to them, it is not accompanied by 
a lattice deformation.0 The characteristic wave number of the CDW component is 2kp, where kp is a Fermi wave 
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This density component is usually taken into account when the system contains an impurity.Lra The latter pins 

S CDW, resulting in the electron density oscillations around the impurity, and in the dc conductivity suppression. It is 
obvious, that since short-range electron correlations contribute to the density, they must also contribute to transport 
■ even if any impurities are absent. In the present paper we investigate this contribution. 

The first problem we meet here is that a conventional bosonized density operator, describing CDW, does not 
conserve the number of particles in the system, and, consequently, fails to describe correctly Id electron transport. 
We propose an expression for the density operator that is free from the mentioned shortcoming and calculate CDW 
contribution to the density response function. The imaginary part of the CDW susceptibility exists in the narrow 
band of wave numbers near 2kp, whereas the smooth component produces a <5-peak in the long-wave region. Both 
CDW and smooth terms give, generally speaking, comparable contribution to the dissipative conductivity, which can 
be determined via the absorbed power. 

In order to estimate CDW contribution to dissipative conductance, the following scheme of measurement was 
proposed. We consider the power dissipated in a Id electron system, subjected to the local external potential that has 
a sufficiently large Fourier-component at 2k p. Short-scale field of this kind can be set up in reality, using near-field 
probing microscopy technique. Our calculations lead to the unexpected conclusion, that CDW contribution to the 
conductivity dominates in the low-frequency limit over that of the commonly discussed smooth component of the 
density. 
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II. DENSITY OPERATOR IN THE LUTTINGER MODEL 



In this section we show that a standard form of the density operator, describing CDW, is not consistent with the 
particle number conservation or, equivalently, with the total electroneutrality requirement. Hence, it is to be altered. 
We aim to obtain an adequate expression for the Id electron density operator. 

Particle number conservation requires that 



+00 



dx (p(x)) = 



where (p(x)) is a density fluctuation at the position x. This requirement can also be formulated for the susceptibility 
Xquj of the system to external potential. The susceptibility Xgw, defined by the linear relation between the Fourier- 
transform of the density {p) gu) and the external potential ip qu , via 
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must have the following limit when w^O 
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In the spinless Luttinger model, the electron density operator p(x) is usually taken aa 



1 kp 

p(x) = d x <p H cos(2fcF2i — 2<j>) 
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where cf>(x) is a bosonic phase. 



-1 


S' 


r=+l. 














\ K k 




(a) 



(b) 



FIG. 1. (a) The Luttinger model spectrum. Arrows show electron transitions, resulting in the smooth (q ~ 0) and CDW 
(q ~ 2/cf) components of electron density, (b) Electron-hole pair spectrum. Parabolas restrict the darked area of all possible 
energies of electron-hole pairs. The low-energy pair excitations exist only in two regions, shown in black, where q ~ and 
q ~ 2k F . 



Two components in the electron density operator correspond to two kinds of the low-energy pair excitations in a 
Id system, shown in Fig.^(a). The first component describes long- wave fluctuations, which are smooth on the kp 
scale and small as compared to background density kp/ir. This part of the electron density operator represents the 
sum of the densities of right- and left-moving r-fermions (denoted by r = ±1).0 It is essential, that the corresponding- 
pair excitations come about separately within each branch of the Luttinger liquid spectrum, and transfer momentum 
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q <C hp. The second (CDW) component of the electron density operator is due to the interference of the left- and 
right-moving electrons. This term describes excitations, transferring momentum ~ 2kp. 

It is easily seen that the form of the smooth density component p sm (x) obeys the particle number conservation 
requirement. The integral of p sm (x) is zero in consequence of boundary conditions: 



+00 ■[ 
dxp sm (x) = <f>(x) 
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The corresponding susceptibility Xam(q,u) indeed goes to zero, when q — > 0, as can be seen from the expression for 
Xsm(?,w) given in Ref. |] (see also Eq. (^) below). However, the CDW part of the electron density operator does not 
exhibit such a property! The integral of p c dw{x) is not zero: 
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dx pcdw(x) = — / dx cos(2kpx — 2(f)) ^ 
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This implies that the integral of the average density (p c d-w{x)) over 2, generally speaking, is not zero too. It is for this 
reason that the susceptibility Xcdw(<7i w ) calculated with the density operator (|l|) does not go to zero when q — > 0, as 



we show in Section HI 



Consequently, the above expression for p c d w and the corresponding density response function x c dw contradict to 
particle number conservation and do not correctly describe Id electron transfer. 

In order to understand the reasons of this failure, we ought to turn to the derivation of Eq. (|l|) for the electron 
density operator in the Luttinger model. 

The field operator ^ e \(x) of electron system is known to be non- locally connected with that of r-fermion one ^(a^.cl 
However, in the low-energy approximation it is usually taken as 

* c i(x) ~ V+(x) + . 

This local expression for a field operator leads to the electron density operator of Eq. (|]) , which does not conserve 
the number of particles. 

The electroneutrality is violated because the consistent low-energy expansion of ^oi^) has not been done. In order 
to find the proper form of the operator, one can either use the correct nonlocal expression for ty c \ or try to find an 
adequate correction to the density operator of Eq. ([!]). We use the second option. Keeping in mind that the density 
operator is a model one, we bring it to the exact differential form by adding the lower order term in the low-energy 
expansion: 

Pcdw(z) = — ( 1 - J cos(2fc F x - 2$) 
n \ k F J 

(2) 

= —d x s\x\{2kpx — 2(f) 

Z7T 

This form of the CDW operator guarantees the number of particles to be conserved. The integral of corrected p c( jw 
is zero in consequence of boundary conditions, and the susceptibility has a correct long-wave limit: Xcdw(9,w) ~ q 2 
when q — > 0. I_ . 

Another derivation of the density operator, conserving the number of particles, was given by Haldane.EHl He suggested 
to take a function (f>(x) that increases by it on each encountered electron. Then the electron density operator is 
presented as 




P(*) = -5Z *+> ^7— ■ (3) 



If we associate <f>(x) with kpx — <fi(x), then the first term of Eq. (g) gives the smooth part of the density in Eq. (Q). 
The first harmonic (m = ±1) of Eq. (|3|) coincides with Eq. (||) in form, but differs by a factor of two. 

In Section [H we show that the density response, calculated with operator of Eq. (g), gives the correct transition to 
the limiting case of noninteracting electrons. By contrast, Eq. (||) and the conventional expression (|l|) for the CDW 
operator fail to give this transition. 

Thus, the bosonized density operator (|^) of CDW component is consistent with particle conservation and can be 
used to describe correctly the density response of interacting electrons. 
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III. CDW CONTRIBUTION TO TRANSPORT 



In this section the contribution of the CDW and smooth components of the density to the linear response function 
and to the dissipative conductance is explored. To be specific, we consider the system, illustrated by Fig. ||, that 
comprises the gated quantum wire and the conducting tip. The tip allows one to act locally on the electrons in the 
wire by ac-potential. Such a configuration is very interesting in view of the rapid progress of near-field microscopy 
technique- This system is close to the metal-carbon nanotube-tip structure considered recently by Y. Xue and 
S. DattaB. 

We assume that the following hierarchy of the scales takes place: 

a < d < kp 1 < I , 

where a is the diameter of the quantum wire, d is the distance between the wire and the gate, and / is the distance 
between the tip and the wire. The first pair of inequalities allows one to consider electrons in the wire as a Id system 
with short-range interaction. The right inequality means that the influence of the tip is much weaker than that of the 
gate. Thus the presence of the tip does not affect electron-electron interaction and the ground state of Id electrons 
is uniform. 
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Quantum wire 
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Gate 



FIG. 2. A view of the model Id system with short-range interaction. The tip produces a local potential disturbance. 



When ac- voltage is applied to the tip, the electrons in the wire are subjected to the electric potential <p(x), which 
is localized in space: (p(x — > ±00) — > 0. If tp(x, t) = ip(x) cos cut, the power dissipated in the system is shown in 
Appendix A to be 

P = -e 2 u> / +0 °^|V(^) • (4) 
Jo 27r 

The conductance of a mesoscopic structure is determined usually through the current and the voltage applied across 
the leads. Since no current-carrying leads are taken into our consideration, it is the dissipated power that determines 
the conductivity via a = P/V 2 , where V is an external potential amplitude in the wireu 

Let us study the contribution of both electron density components to the susceptibility and conductivity. The 
density response function is found via Kubo formula with the density operator p = p sm + p c d w , Pcdw being defined in 
Eq. (§. 



A. Long-wave density fluctuations 

The density response function, corresponding to p sm , is well investigated.0 Nevertheless it is given below to demon- 
strate that the total electroneutrality is kept, and to compare the dynamics of smooth density response with that of 
CDW: 

X S m(x,t) = l-6(t)d x (6(vt + x)-6(vt-x)) , (5) 
h 
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where g is the interaction parameter and v is the velocity of boson excitations .0 It is obvious that the integral of 
Xsm(%,t) over x is zero at every moment of time. The susceptibility 

Xsm {q,u}) = -- — — — 5-2 ■ W 

h (lo + iQ) 2 — q z v z 

goes to zero, when q — > 0. This means that the number of particles is conserved. 

Eq. (||) shows that the evolution of the electron density fluctuation, generated by the local disturbance of the form 
ip(x,t) ~ S(x)5(t), is presented by two wave fronts propagating in opposite directions with a constant amplitude and 
the velocity v, as shown in Fig. 0. These waves transfer charge and produce current. 
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FIG. 3. The response function Xsm (a;, t) for a given t. 



The fact that electron current depends on the position x along the Id conductor should not be perceived as a 
confusing one. It does not contradict to the general requirement of the electrodynamics according to which the total 
current / must be conserved, i.e. be independent of x. The current conservation is obeyed if we take into account 
the displacement current. It was shown in Ref. ^ that the displacement current between the Id conductor and the 
gate ensures the total current conservation. 

The dissipated power due to the smooth density component is 



e 2 w 2 



2h v< 



(7) 



where (p u / v is the Fourier-component of tp(x) at q = uj/i 



B. 2k p density fluctuations 

Now we turn to CDW contribution to electron transport. The response function, calculated with the density 
operator is as follows: 



Xcdw{x,t) 



where a = kp 1 , and 



9{t) 3 2 
2irh dx 2 



a 2 + (vt + x) 2 a 2 + (vt — x) 2 
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in (Trg8 a (v 2 t 2 - x 2 )) cos(2/c^x) 



(8) 



n9 a (v t — x ) = arctan 



vt + x 



arctan 



vt — x 
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9 a (z) is a smoothed step function that goes to the Heaviside -function, when a — > 0. It describes the smoothed 
fronts, moving with the velocity v. 

The integral of Xcdw(x, t) over x is zero for any time t because of the presence of the second derivative. The response 
function, calculated with the conventional expression (Q) for p c d w , differs from Eq. (|^) just in the absence of the second 
derivative with respect to x. The integral of this response function is obviously non-zero. 

The response function Xcdw (x, t) represents the rapid oscillations modulated in the amplitude by the envelope that 
has the form of the moving wave, as illustrated in Fig. ^. The wave front propagates with the velocity v, and the 
amplitude of the wave diminishes with time as t~ 2g . 
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FIG. 4. The CDW response function for a given t. 



The CDW susceptibility x c d w (<Z,^) is 

Xcdw(g,w) = - 



ft w 45+ir 2 ( 5 )sin(7rg) \k F 
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where x'i (<3S w ) is a rea l function without any singularities. The function x'i (<Zi w) is not important if we are interested in 
the dissipative conductance. Hence, the imaginary part of the CDW susceptibility is expressed in terms of elementary 
functions. Inside the interval \\q\ — 2kp\ < \ui\/v, 



Xcdw (9, w) = 



1 sign u 



vkp 



kF 
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(9) 



and outside this interval, Xcdw(9i a; ) = 0- The details of the derivation of Eq. (^) are given in Appendix B. Equa- 
tion (^) is valid for any interaction parameter g. It is worth noting that, when g = 1, Eq. (^) exactly reproduces 
the corresponding part of the Lindhard formula for free electrons. Incidentally, the susceptibility of free electrons, 
calculated with the density operator of Eq. (||), differs from that given by the Lindhard formula by a factor. 

The dependence of the imaginary part of the CDW susceptibility on the wave number q is shown in Fig. [5] for several 
values of interaction parameter. It is seen that, in contrast to the nonintcracting case, the inclusion of electron-electron 
interaction leads to the qualitative change in the shape of x"dw(9j w ) an d to the increase of Xcdw(9> w ) with interaction 
strength. 

The real part of the susceptibility, X c dw(9i w )' ^ s non-zero for all lu and q. When q goes to zero, X c dw(9' w ) diminishes 
as q 2 , in accordance with the exact differential form of the expression (||) for the CDW density operator. Had we not 
changed /O c dw to the exact differential form, we would have gained a violation of the total electroneutrality. 

The qualitative effect of the electron-electron interaction on the CDW susceptibility (on both the real and the 
imaginary parts) consists in the singularity that appears at 2kp ± uj/v. This singularity is obviously a dynamic 
analogue of Kohn's anomaly. 
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FIG. 5. The imaginary part of the CDW susceptibility as a function of the wave number for the different values of the 
interaction parameter and htv/ep — 0.2. 



The dissipated power due to CDW is 



P< cdw — 



h A9T{g)T(g+\) \vk 



2</ 



kp \if2la. 



(10) 



where <p2fe F is the Fourier-component of the external potential at q — 2kp. The dependence of P c d w on the interaction 
parameter g is shown in Fig. g for two frequencies. P c dw is seen to have a pronounced peak at g ~ ln _1 (eF/ n ^)j 
where £f is the Fermi energy. 
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FIG. 6. The dissipated power P ct j w as a function of the interaction parameter g. E2k F is the Fourier-component of the 
external electric field at q — 2fc_F- 



Let us compare the dissipated power for the CDW and conventional conductivity mechanisms. Equations (0) and 
( |io| ) show that 



P< 



cdw 
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Since g < 1, the CDW response is seen to dominate in the low-frequency limit. 



IV. CONCLUSION 



In the present paper we have investigated the effect of short-range electron correlations on the dynamic conductivity 
of a Id conductor within the frame of the Luttinger model. Short-range electron correlations are described by the 
CDW component of the bosonized density operator. We have found that a conventional form of this operator is 
not adequate to describe electron transfer. Namely, it does not conserve the number of particles in the system. Wc 
relate this inconsistency to the fact, that the low-energy expansion of this operator is not properly performed in 
the Luttinger model. We have proposed a model operator of electron density that does not contradict to the total 
electroneutrality requirement and gives the correct transition to the case of noninteracting electrons. 

The CDW density response function represents the rapid oscillations, modulated by the envelope function, which 
has the form of a moving wave. The CDW density component contributes significantly to the susceptibility of the Id 
conductor without any impurities. The imaginary part of the CDW susceptibility is non-zero in the band of width 
w/v around 2kp- Within this band the dissipative susceptibility depends strongly on the interaction parameter g. 

We have proposed a scheme of experiment in which the CDW mechanism of electron transport dominates over 
the commonly discussed mechanism due to long-wave fluctuations. The measurement can be realized with the use 
of the local excitation by the tip with ac-potential. In the low-frequency regime the power dissipated in the system 
is determined predominantly by the CDW transport mechanism. The power drastically depends on interaction, 
providing an opportunity to study interaction effects in quantum wires. 
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APPENDIX A: THE POWER DISSIPATED IN A GATED WIRE 

The goal of this section is to obtain Eq. (|J) for the power dissipated in the system, shown in Fig. ^. 

When the voltage is applied to the tip, it produces electric potential (p ex t(r,t). The external field induces charges 
in the wire and in the gate, so that the total electric potential is the sum of the external potential and the induced 
one y>ind(r, t). The induced potential is determined by a linear equation 

L (^md = p(r, t) , 

where L is the Laplace, with appropriate boundary conditions at the gate. It coincides with the equation that governs 
the electron-electron interaction potential V^(r): 

LV(r-r') = -<J(r-r') . 

Consequently, the interaction potential is a Green's function of this equation with boundary conditions imposed by 
the gate. So, tp^ is connected with V via 



Wnd(r) =-J dv' V(r - r')p{r') . 



The effective Id potential in the wire <p(x, t) is connected with 3d potential <p(r, t) through the following integral 
relatiorO 

<fi(x,t) = f d 2 r±<p(r,t)g(r_ L ) , 

where g(r±) is a distribution function of electron density along the radial coordinate in the wire. 
Hence, the induced potential in the wire satisfies the equation 

Pindfow) = ~U(q)p{q,uj) , 
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where U{q) is a Id effective interaction potential with account of the gate. Thus the total potential and the external 
one are connected as follows 



<Aot(g,w) = <A=xt(g,w)[i - U(q)x(q,v)] ■ (A.l) 

The dissipated power is given by 



+00 



P= dx E tot (x,t)j(x,t) , 



—00 



where j(x,t) is electron current and E to t(x,t) is the total electric field. Using Eq. (A.l) for the total potential and 
the continuity equation, one obtains the formula (Q) for the power dissipated in the wire. 

The presence of the displacement current in the gap between the wire and the gate implies that electron current 
flows inside the gate. This current is of the same order of magnitude as the electron current in the wire. However, its 
contribution to the dissipated power is negligible because the conductivity of the gate is extremely high. 

APPENDIX B: FOURIER-TRANSFORM OF THE CDW RESPONSE FUNCTION 

CDW linear density response is 

y(x t) m ^ ( * 29 cos(2k F x) 
Xy , } 8ir 2 hdx 2 \[(a + i(vt + x))(a + i(vt - x))]a 

We reduce the calculation of the Fourier-transform to the calculation of one integral: 

+ 00 +OC 



X (q,uj)= J dte^ j dxe^xOM) 
-00 

q {3{q + 2k F ,Lu)+3(q~2k Fl Lu)-3*{q + 2k Fl -Lu)-3*(q-2k F ,-cu)), 
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where the integral to calculate is 

+00 +OC 

3(q,u>)= J dte wt j dxe tqx (f3 + ix)' 9 (P-ixy 9 ; = a + ivt; a, t, g > 0. (B.2) 
-00 

First, we calculate the integral with respect to x: 

+00 

1= J dxe lqx (p + ixy 9 ([3-ix)- 9 . 

— OO 

The integral is an even function of q; first we take q > 0, and then q — -> \q\. The key point is to use two representations 
for a MacDonald functioriij: 

+00 

K v (z) = &Lr\u+^\ J dt (t 2 + z 2 )~ u ~host, Rez/>-~, \argz\ < ~, 



+°° 1/— 1 

K M = \J\ r_1 z " e ~ Z J dte- zt f-i + 2 ,Rez>0, Re^> ~ . 



Then the integral / can be transformed in the following way 
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Taking this form of I into account, let us now calculate the integral over t in Eq. (|B 

+00 
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lu — qv(z + 1) + iO 



- p -q(z+l)a 



it5(lo - qv(z + 1)) + iP 



1 



a; — qv{z + 1) 



A real part of the integral He3(q, lo) may be expressed via elementary functions. (5-function contributes, when its 
argument zq > 0: 

id 

LO - qV(Z + 1) = 0, Zq = 1, 

qv 



that is when q < — , lo > 



4-7T 2 

Re%u;) = — - — r?9 



q-» 1 / e~ qa( - z+1 '>z 9 - 1 (2 + zy- 1 S{uj-qv(z + l))dz 
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Hence, the imaginary part of the CDW susceptibility, determined by Eq. (B.l), is 
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,Qr<~, 

V 



where Q r = \q — r2kp\. 

Bosonization requires Two -C Ep. Therefore 



Xcdw(?» w ) = 



1 1 

T^A9T 2 {g) 



LO 
vkF 



\q\_2kp 
kp 



-1 s-i 



M -2fc F | < - . 

v 



This formula is valid for any g > 0, and when w > 0. When to < 0, x'^QS^) — ~x"(q, M 
The real part of the susceptibilty is 



Xcdwfew) 



1 



1 



47h; A9T 2 {g) sin(7rg) \kp 
where Xiili^) is a rea l function without any singularities 
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